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We report on the analytical calculation of NNNLO {0{a^)) conversion factor between the MS quark mass and 
the one defined in the so-called "Regularization Invariant" scheme. The NNNLO contribution in the conversion 
factor turns out to be relatively large and comparable to the known NNLO term. 



1. Introduction 

Quark masses are fundamental parameters of 
the QCD Lagrangian. Nevertheless, their relation 
to measurable physical quantities is not direct: 
the masses depend on the renormalization scheme 
and, within a given one, on the renormalization 
scale fi. 

In the realm of pQCD the most often used defi- 
nition is based on the MS-schemeJ^]j^ which leads 
to the so-called short-distance MS mass. Such 
a definition is of great convenience for dealing 
with mass-dependent inclusive physical observ- 
ables dominated by short distances (for a review 
see Q). Unfortunately, as their mass dependence 
is relatively weak the predictions are usually dif- 
ficult to use for getting a precise information on 
quark masses. 

To determine the absolute values of quark 
masses, one necessarily has to rely on the methods 
which incorporate the features of nonperturbative 
QCD. So far, the only two methods which are 
based on QCD from the first principles arc QCD 
sum rules and lattice QCD (for recent discussions 
see e.g. P-|^,p|-p^). Rather accurate determina- 
tions of the ratios of various quarks masses can be 
obtained within Chiral Perturbation Theory . 

Lattice QCD provides a direct way to deter- 
mine a quark mass from first principles. Un- 
like QCD sum rules it does not require model 
assumptions. It is possible to carry out the 
systematic improvement of lattice QCD so that 
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all the discretization errors proportional to the 
lattice spacing are eliminated (a comprehensible 
review is given in [p^ ). The resulting quark 
mass is the (short distance) bare lattice quark 
mass. The matching of the lattice quark masses 
to those defined in a continuum perturbative 
scheme requires the calculation of the correspond- 
ing multiplicative renormalization constants. In 
the RI scheme [|l^ the renormalization conditions 
are applied to amputated Green functions in Lan- 
dau gauge, setting them equal to their tree-level 
values. This allows the non-perturbative calcu- 
lation of the renormalization constants. An al- 
ternative to the RI approach is the Schroedinger 
functional scheme (SF) which was used in 

Once the RI quark masses are determined from 
lattice calculations they should be related to the 
MS mass by a corresponding conversion factor. 
By necessity this factor can be defined and, hence, 
computed only perturbatively. The conversion 
factor is now known at next-to-next-to-leading or- 
der (NNLO) from Ref. |l9|. The NNLO contribu- 
tion happens to be numerically significant. This 
makes mandatory to know the NNNLO 0{al) 
term in the conversion factor. 

In this work we report on the calculation of this 
term. It turns out that the size of the newly com- 
puted term is comparable to the previous one at 
a renormalization scale of 2 GeV — the typical 
scale currently used in the lattice calculations of 
the light quark masses. This means that pertur- 
bation theory can not be used for a precise con- 
version of the presently available RI quark masses 
to the MS ones. A simple analysis shows that the 



convergency gets better if the scale is increased 
to, say, 3 GeV. Thus, once the lattice calcula- 
tions produce the RI quark masses at this scale 
our formulas will allow accurate conversion to the 
MS masses at the same scale. 

2. Scheme dependence of quark mass 

We start by considering the bare quark propa- 
gator (for simplicity we stick to the Landau gauge 
and, thus, do not explicitly display the gauge de- 
pendence) 



Fo(g,a°,TOo) = (mo - ^ - So) \ 



(1) 



with the quark mass operator So being conve- 
niently decomposed into Lorentz invariant struc- 
tures according to 



So =^S?. + TOoS!^. 



(2) 



Here mo and Tpo is the bare quark mass and field 

respectively: a. . ^ = ^ and 5 is the bare 
QCD gauge coupling. To be precise we assume 
that (|l|) is dimensionally regulated by going to 
non-integer values of the space-time dimension 
D = 4 - 2e [|0|j2|]. The MS renormahzed coun- 
terpart of the Green function (M reads 



F(g,as,m,/i) = (m - ^ - S) ^ 
^Z-'Fo{q,alm,o)\ 



(3) 



— 1/2 

where the renormahzed quark field — Z2 wo 
and the 'Hooft mass parameter /i is a scale at 
which the renormahzed quark mass is defined. 
The renormalization constants Z2, Za and Zm are 
series of the generic form ( ? = 2, a or m ) 

.si , N , /n \ i calculat 

^? = 1 + E^'^^' zW=^Z^)piy.(4) theMS 



i>0 



j>i 



Now let us consider the quark propagator renor- 
mahzed according to a different subtraction pro- 
cedure. Marking with a prime parameters of the 
second scheme one can write 



F{q,a'^,m',fi) = - 
= {Z!,)-^Fo{q,almo)\ 



1 



(5) 



where without essential loss of generality we have 
set /i' — fi. The finiteness of the renormahzed 
fields and parameters in both schemes implies 
that, within the framework of perturbation the- 
ory, the relation between them can be uniquely 
described as follows 



m 



Cn 



(6) 



with the "conversion functions" being themselves 
finite series in a^, i.e. 



Cv = 1 



i>0 



(7) 



for ? = m or 2. 

Note that in general the coefficients C7 may 
depend on the ratio m! j [i. If such a depen- 
dence is absent then the corresponding subtrac- 
tion scheme is referred to as a "mass indepen- 
dent" one. In what follows we limit ourselves 
to considering the latter case. In addition, be- 
ing only interested in conversion functions C2 
and Cm, we will assume that the function Cq 
has already been determined and, thus, will deal 
with the following representation of C2 and Cm in 
terms of the MS coupling constant as ( ? = 2,m 
) 



C7 EE 1 



i>0 



(8) 



From eqs. ((|i|)) it is easy to see that 

C2-(l + Sv) = l + SV, (9) 
C2 • Cm • (1 - Ss) = 1 - S'5. (10) 

The renormalization conditions for the non-MS 
scheme should then be used to provide the nec- 
essary information about the right hand side to 
calculate the conversion factors C,„ and C2 given 
the MS renormahzed Sy and S5. 

An example of a mass independent MOM 
scheme has recently been suggested in un- 
der the name of RI ( "Regularization Invariant" ) 
scheme. The corresponding versions of the 
general conversions formulas read ||2^ ( I = 



C^' 



asy(^) 



dl 



(11) 
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3. Three Loop Conversion functions 



We have analytically computed the functions 
Tiy and in the massless limit to order a^. 
The calculation has been done with intensive use 
of computer algebra programs. In particular, we 
have used QGRAF f or the generation of dia- 
grams, and MINCER [^4[ for their evaluation. 

Our results for the conversion functions read 
(in the Landau gauge and in terms of fractions 
and Riemann's Zeta function and as functions of 
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At a scale — 2 GeV and Uf — 4, the numer- 
ical contributions of various orders are as follows 
(for simplicity we took as(2 GeV)/7r = .1) 



C5 = 1. - 0.1333 - 0.0754 - 0.0495 



and 



1. + 0.- 0.00477- 0.00508. 



(13) 



(14) 



One observes that the sizes of the NNLO 
and NNNLO contributions to C^^^ , at this scale 



amount to about 7.5% and 5% respectively. This 
means that perturbation theory can not be used 
for a precise conversion of the RI quark masses to 
the MS ones at the renormalization scale /i = 2 
GeV. The convergency can be improved if one 
icreases fj, till, say, 3 GeV. Indeed, with this 
choice of /i the standard three-loop evolution 
gives asiS GeV) = 0.262 and eqs. ^j^) trans- 
form to 



C^/ = 1. - 0.111 - 0.0526 - 0.0289 
and 

C^' = 1. - 0.00333 - 0.00296. 



(15) 



(16) 



To conclude: the newly computed NNNLO cor- 
rections are numerically significant and should 
be taken into account when transforming the 
RI quark masses to the MS ones. 
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